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ABSTRACT 



The supernova SN 2001du was discovered in the galaxy NGC1365 at a distance of 
19 ± 2 Mpc, and is a core-collapse event of type II-P. Images of this galaxy, of moder- 
ate depth, have been taken with the Hubble Space Telescope approximately 6.6 years 
before discovery and include the supernova position on the WFPC2 field of view. We 
have observed the supernova with the WFPC2 to allow accurate differential astrome- 
try of SN 2001du on the pre-explosion frames. As a core-collapse event it is expected 
that the progenitor was a massive, luminous star. There is a marginal detection (3cr) of 
a source close to the supernova position on the prediscovery T^— band frame, however 
it is not precisely coincident and we do not believe it to be a robust detection of a 
point source. We conclude that there is no stellar progenitor at the supernova position 
and derive sensitivity limits of the prediscovery images which provide an upper mass 
limit for the progenitor star. We estimate that the progenitor had a mass of less than 
15Mq. We revisit two other nearby SNe II-P which have high quality pre-explosion 
images, and refine the upper mass limits for the progenitor stars. Using a new distance 
determination for SN 1999gi from the expanding photosphere method, we revise the 
upper mass hmit to 12M0. We present new HST images of the site of SN 1999em, 
which validate the use of lower spatial resolution ground-based images in the pro- 
genitor studies and use a new Cepheid distance to the galaxy to measure an upper 
mass limit of 15Mq for that progenitor. Finally we compile all the direct information 
available for the progenitors of eight nearby core-collapse supernovae and compare 
their mass estimates. These are compared with the latest stellar evolutionary mod- 
els of pre-supernova evolution which have attempted to relate metallicity and mass 
to the supernovae type. Although this is statistically limited at present, reasonable 
agreement is already found for the lower mass events (generally the II-P), but some 
discrepancies appear at higher masses. 

Key words: galaxies: individual (NGC1365) - stars: evolution - supernovae: general 
- supernovae: individual (SN 2001du, SN 1999em, SN 1999gi) 



1 INTRODUCTION 



The end point of the lives of all stars with initial mass greater 
than approximately SM© is thought to be a core-collapse su- 
pernova. Supernovae (SNe) explosions which arise from mas- 
sive stars have been a major driver in the chemical evolution 
of the Universe, are important in shaping the interstellar 
medium dynamics in gas-rich galaxies and are fundamental 
in the study of the origin of the chemical elements in the 
Universe. The exact evolutionary states of stars before they 



undergo core-collapse is of great interest in the fields of mas- 
sive star evolution and the physics of supernovae evolution. 

Supernovae are classified according to the atomic lines 
observed in the early-time optical spectra. Briefly, the pres- 
ence of Hi optical lines indicate a type II classification 
(SNe II), while those that show no hydrogen at all are 
type I (Filippenko 1997). The SNe la are thought to arise 
through thermonuclear explosions in white dwarf binary sys- 
tems (Branch et al. 1995). Given the occurrence of these 
supernovae in both old and intermediate age stellar popu- 
lations (e.g. in galaxies of type SO and earlier), the white 
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dwarf binary scenario is consistent with the progenitors be- 
ing low— intermediate mass stars. However very little other 
details are known about the stellar systems that produce 
SNe la. All other supernovae including the type Ib/Ic (SNe 
Ib/Ic) and all the subtypes of SNo II arc thought to be due 
to core-collapse during the deaths of massive stars. The sub- 
type SNe II-P show prominent, broad H I lines with P-Cygni 
type profiles in their optical spectra, indicating that the 
progenitor retained a substantial hydrogen envelope prior 
to explosion e.g. see Woosley & Weaver (1986). They also 
display a plateau phcise in their lightcurve during which the 
supernovae remain at approximately constant brightness for 
a period of between 30 and 100 days, which is also inter- 
preted as evidence of an extended hydrogen atmosphere in 
the progenitor star, such as red supergiant (Chevalier 1976, 
Woosley & Weaver 1986). This plateau phase is most obvi- 
ous in the l^— band, but can less pronounced or non-existant 
at other visual and NIR wavelengths. The type-IIn super- 
novae show strong hydrogen features which are in emission 
with broad wings and very narrow cores. These features are 
thought to arise from interaction with dense circumstellar 
material which may have been ejected from the progenitor 
star through stellar winds or more dramatic mass-loss events 
e.g. see Fassia et al. (2001), Pastorello et al. (2002), Lentz 
et al. (2001) and references therein. The SNe Ib/Ic also do 
not show any obvious signs of hydrogen in their spectra, 
although SNe lb display pronounced He I absorption. How- 
ever recently Branch (2003) has suggested that at least some 
lb events do show some amounts of hydrogen, and some Ic 
events do show helium. It is likely there is a continuum of 
possibilities rather than the well defined sub-classes tradi- 
tionally employed. 

There is very strong evidence that the SNe II and Ib/Ic 
are associated with the deaths of massive stars as they are 
never seen in elliptical galaxies, only rarely in SO types (van 
den Bergh et al. 2002) and they often appear to be associ- 
ated with sites of recent starformation such as Hll regions 
and OB associations in spiral and irregular galaxies (Van 
Dyk et al. 1996). However the desired direct evidence of de- 
tecting stars before they exploded has only been possible 
for a small number of objects. The two most unambiguous 
detections of a progenitor object of confirmed supernovae 
explosions are for SNe 1987A and 1993 J. The progenitor to 
SN 1987A in the LMC (White & Mahn 1987) was the blue 
supergiant Sk— 69°202 of spectral type B3Ia (Walborn et al. 
1989). The closest supernova to the Milky Way since then 
was SN 1993J in M81 (3.63 Mpc), which was a type lib event 
and ground based U BVRI photometry of the SN site be- 
fore explosion was presented by Aldering et al. (1994). The 
photometry of the progenitor candidate was best fit with a 
composite spectral energy distribution of a KO la star and 
some excess UB band flux either from unresolved OB as- 
sociation contamination or a hot companion. Neither the 
progenitor of SN 1987A nor that of SN 1993J is consistent 
with the canonical stellar evolution picture, where core car- 
bon burning finishes and core-collapse occurs relatively soon 
afterwards (~ 10^ — 10* yrs) while the massive star is an M- 
supergiant. Claims have been made on the detection of the 
progenitors of three other supernovae events SNe 1961V, 
1978K and 1997bs, however for reasons we discuss in Sec- 
tion 5.3 it is uncertain whether these were actual supernovae 
explosions. 



The only way to further add to our knowledge of the 
core-collapse process and its relation to massive stellar evo- 
lution is to identify more progenitor objects. There is now a 
wealth of data in the Hubble Space Telescope (HST) archive, 
as well as high quality ground-based images from various 
weU maintained archive facUities (e.g. ESO, ING, CFHT). 
Thus there are reasonable chances of having available an 
image of a pre-supernova site which has suflScient resolution 
to identify individual massive stars. Clearly there is a dis- 
tance limit within which images are useful from depth and 
resolution constraints. This is roughly 10 Mpc for a deep, 
good quality (natural seeing of typically 0.7") image from a 
ground-based 2-8m telescope, and roughly 20 Mpc for typical 
exposure length HST frames. Recently Smartt et al. (2001a; 
2002a; 2002b) have used both HST and ground based im- 
ages of pre-explosion supernovae sites to try to detect pro- 
genitor stars of SNe II-P 1999gi and 1999em and the type 
Ic SN 2002ap. These nearby events have high quality pro- 
discovery archive images available, but no progenitor star 
is actually identified in any of these cases. These papers 
discuss the luminosity limits which can be placed on the 
exploding star from the sensitivity measurements of the im- 
ages, and quite restrictive mass-limits can be determined 
(see Section 5.3 for a full comparison and discussion). In a 
similar vein. Van Dyk et al. (2003a) have attempted to iden- 
tify the progenitors of 16 type II and type Ib/s SNe on HST 
prediscovery images. However 9 of the objects in this lengthy 
paper are further than 'N^40Mpc, with another reddened by 
some 5 — 6 magnitudes of visual extinction, hence there is 
little or no chance of identifying a progenitor star for these 
events. The other six are subject to the detailed HST follow- 
up programme of Smartt et al. (GO 9353) which will allow 
more accurate astromctry and possible source detection. The 
Van Dyk et al. study has attempted absolute astrometry of 
WFPC2 archive images as a method of determining the su- 
pernovae positions on the prediscovery frames. This method 
has not proved accurate enough to position the events to 
the precision necessary to assess if a progenitor star can be 
detected. The error boxes the define are generally too large 
to provide useful constraints. They do sugges that there are 
possible progenitor candidates for four supernovae, although 
these will require more detailed follow-up to confirm or dis- 
count the claims. 

In this paper we present a search for a progenitor of 
the core-collapse supernova SN 2001du. This weis a type 
II-P event in the galaxy NGC1365 at a distance of (m — 
M) = 31.39 ± 0.26 (Ferrarese et al 2000; corresponding to 
19±2 Mpc) and the supernova position was observed by HST 
in 1995. The supernova was discovered by Evans (2001) on 
2001 August 24.7 UT at a visual magnitude of about 14.0. 
It was classified by Smartt et al. (2001c) on 2001 September 
2.05 UT, and Wang et al. (2001) on August 30, as a type 
II-P supernova due to the presence of strong Hi P-Cygni 
profiles in the H/3 and Ha transitions. Wang et al. also ob- 
tained spectropolarimetry of SN 2001du and reported the 
supernova to be polarized at a level of 0.29% with a flat 
spectrum characteristic of electron scattering, indicative of 
a fairly spherical explosion. This supernova has very similar 
optical spectra to the other type II-P events SNe 1999em 
and 1999gi for which we determined upper mass limits for 
the progenitor star. Images of the SN 2001du host galaxy 
are available in the HST archive allowing another detailed 
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search for a precursor object. In this paper we show that 
there is a marginal detection of a faint source close to the 
supernova position. But it is neither a robust detection of 
a single point source nor precisely coincident with the su- 
pernova. Wc argue that there is no detection of a progeni- 
tor, and hence determine mass and luminosity limits of the 
exploding star and discuss the implications of this finding 
alongside all the available, direct information we have on 
progenitors of core-collapsc supernovae. 

2 OBSERVATIONAL DATA AND ANALYSIS 

2.1 Astrometry and detection limits of WFPC2 

images 

The galaxy NGC1365 was observed on four epochs with the 
WFPC2 on board HST before the supernova 2001du was 
discovered. However on only one of these occasions does the 
supernova position actually fall on any of the WFPC2 chips. 
The galaxy was a target for the HST key project on the ex- 
tragalactic distance scale, which was designed to accurately 
resolve and photometrically measure Cepheid variables in 
this spiral. Unfortunately these deep exposures (totalling 
18.5hrs in F555W and 6hrs in F814W) do not cover the 
position of the progenitor of SN 2001du. The observations 
which do include the supernova position were taken in Cy- 
cle 4 in G05222 (PI: J. Trauger), and the details arc listed in 
Table 1. These data were located and retrieved from the ST- 
ECF Science archive using the astrovirtel initiative. The 
approximate position of the supernova on these prcdiscovcry 
images can be estimated (Van Dyk et al. 2003a), however 
one is limited by the absolute astrometric resolution of the 
WFPC2 frames (usually of the order 1 — 2") and the error 
on the SN position (~ 0.5 — 1"). As this could amount to 
~30 pixels on the WF chips of WFPC2, it is essential that 
more accurate astrometry is achieved to determine the SN 
position on the prediscovery frames. As discussed in Smartt 
et al. (2001a), by far the best method is to take repeat 
WFPC2 observations of the supernova some months after 
explosion and use the surrounding stars in the post and pre- 
explosion images to determine precise differential astrome- 
try. The crowded host region of SN 1999gi shown by Smartt 
et al. (2001a) and Li et al. (2002) illustrates that an RMS er- 
ror of even 0.2" in positioning the supernova would not have 
allowed two nearby, resolved point sources to be have been 
ruled out as the progenitor objects. Hence re-observing the 
SN at the resolution of HST (with either WFPC2 or ACS) 
is essential to accurately determine the SN position. 

The supernova was observed with WFPC2 on 2001 
November 26 through the program GO9041, and the ex- 
posure times and filter details are listed in Table 1. These 
are total exposure times, and as the SN was still quite 
bright at this epoch, the F336W, F555W and F814W were 
broken down into sub-exposures to allow photometry of a 
non-saturated SN point-sprcad-function (PSF) as well as 
the fainter background stars. Both the F555W and F814W 
datasets are composed of a 40s and 3 x 200s exposures, while 
the F336W is composed of 100s and 460s exposures (with 
CR splits for cosmic ray removal) . Photometry of stars visi- 
ble in both these long and short exposures showed excellent 
alignment to better than 0.1 pixel in each dimension. Hence 



the centroid position of the supernova can be measured in 
each of the short frames, and assumed to be appropriate for 
the longer frames. Seventeen bright (V555 ~ 20.5 — 23.5), sin- 
gle stars were identified in a 500 x 500 pixel region around 
the SN position in the F555W post-explosion 600s expo- 
sure and their positions measured with aperture photometry 
within the IRAF DAOPHOT package. The positions of these 
same 17 staxs were identified in the prediscovery frames and 
a geometric solution mapping these positions to the post- 
explosion frame was calculated. The spatial transformation 
fitted a linear shift, magnification factor and a rotation, giv- 
ing an RMS to the fit of 0.01" in each of x and y. Eleven 
stars were used in the F814W transformation, with similar 
RMS errors. In Fig. 1 wc show the position of the supernova 
in the 3 colour prediscovery images. The total exposure time 
available for all the stacked images at this position is 330s 
in each filter. However this is made up of 3 x 100s and 3x 10s 
individual frames. Given the significant road-out noise of the 
WFPC2 chips (7e~ per frame) compared with the sky back- 
ground in a stellar aperture, the 3xl0s frames were ignored, 
and only the 300s stacked images were used. We note that 
Van Dyk et al. (2003a) quote an exposure time of 100s total 
for these images, and presume this is a typographical error. 

Very close to the position of the supernova on the 
F555W imago there is a possible detection of a faint source 
which is close to the detection limit of the image. Assum- 
ing Poisson statistics, a total CCD readout noise of 12e~ 
(from the three individual stacked frames), and a gain of 
14e~/ADU, this object is detected with a 3. Oct significance 
above the background noise. The magnitude of the object 
estimated by placing a 2-pixel radius aperture at the po- 
sition is V555 — 25.0 ± 0.4. This was obtained using the 
IRAF DAOPHOT package, and applying the appropriate 
charge-transfer-efficiency (Whitmore et al. 1999) and aper- 
ture corrections (Holtzman et al. 1995a). A model point- 
spread-function was calculated using the Tiny Tim software 
(Krist & Hook 1999), and PSF fitting photometry was also 
carried out on the source using the DAOPHOT task ALL- 
STAR. A consistent magnitude of V555 ~ 25.2 ± 0.4 is mea- 
sured. However it is a marginal detection, the object is not 
clearly identified by eye on the frames and is not a con- 
vincing point source. Placing a 2 pixel aperture centred on 
various pixels in the vicinity which are individually more 
than 2a above their neighbours will frequently give a mea- 
sured magnitude of between 25.0—26.2. The source appears 
somewhat extended, and has two peats of intensity sepa- 
rated by roughly two pixels. After PSF subtraction one of 
these pixel peaks always remains, which suggests to us that 
this is not a robust detection of a point source. The source 
is not detected in the F814W filter. Furthermore, the su- 
pernova position is not exactly coincident with the object 
within the estimated errors of our methods. We estimate 
the error radius on the precursor object position to be 0.3 
pixels, which is the maximum difi'erence between positional 
determinations using four methods i.e. intensity weighted 
means in each dimension, a Gaussian fit, a triangular weight- 
ing function, and the PSF fit. However we find the object 
1.03 pixels from the supernova centroid. The error in the 
supernova centroid is 0.14 pixels, and transformation error 
as discussed above is 0.14 pixels RMS, giving a total error 
circle of radius rerr = 0.28 pixels. Hence considering the er- 
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Table 1. Details of the HST WFPC2 archival and GO observational data for SN 2001du. The dataset names come from the the ST-ECF 
Science archive (http://archive.eso.org) and are the associations of shorter, aligned exposures. The location column indicates where 
SN 2001du occurred on the WFPC2 FO\'. The PCI has a T)lato scale of 0.0155" pn- jjixol. WF3 has 0.0996" per pixel. 



Dataset 


Date 




Filter 


Exposure Time 


Location 


U2KV010AB 


1995 


January 15tli 


F160BW 


330s 


WF3 


U2KV010AB 


1995 


January 15th 


F336W 


330s 


WF3 


U2KV010JB 


1995 


January 15th 


F555W 


330s 


WF3 


U2KV010SB 


1995 


January 15th 


F814W 


330s 


WF3 


U6BR030JB 


2001 


November 26th 


F336W 


560s 


PCI 


U6BR030NB 


2001 


November 26th 


F439W 


400s 


PCI 


U6BR0303B 


2001 


November 26th 


F555W 


640s 


PCI 


U6BR030LB 


2001 


November 26th 


F675W 


400s 


PCI 


U6BR030BB 


2001 


November 26th 


F814W 


640s 


PCI 



Figure 1. Panels (a)-(c) show images of the F336W, F555W and F814W filters taken 6.6 yrs before explosion transformed to the pixel 
coordinates of the post-explosion (d). The images (a)-(c) are from the WF3 chip (0.1" pixels), and (d) is better sampled with the 0.05" 

pixels of the PCI. The position of the SN is marked as a cross in the prediscovery images, and the blank triangular regions arc due to 
the edge of the WF3 chip after transformation. The fields of view of these panels are 20", corresponding to a distance of 1.8 kpc at the 
distance of NGC1365 



rors in the analysis, the supernova position is not precisely 
coincident with this possible source (see Fig. 2). 

We have followed the labelling of stars in the region from 
Van Dyk et al. (2003a), who identify three objects within a 
0.9" error circle of their estimated supernova position, but 
the only unambiguous detection we confirm is that of star 
A. The object we marginally detect close to the supernova 
position is the Van Dyk et al. star B, and we keep the nomen- 
clature for clarity. Their star C is also a marginal detection, 
but as it is some way from the supernova position, we do 
not consider it any further. These images are available in 
the HST archive, and we quote our determination of the po- 
sition of SN 2001du (at pixel (106.80, 675.96), with errors 
as listed above) on chip WF3, to allow our results to be in- 
dependently verified if necessary. There doesn't appear to 
be a significant discrepancy between our results and those 
of Van Dyk et al., as the latter have used the HSTphot 
photometry routines of Dolphin (2000a; 2000b) with a 3a 
threshold. We do repeat a 3a detection, but for the reasons 
outlined above do not interpret it as a certain detection of 
a stellar progenitor. Van Dyk et al. measure a magnitude of 
25.02±0.3, which is in good agreement with our PSF mag- 
nitude given the errors, the faintness of the source and the 
different methods employed. 

We determined the detection limits of the data in the 
vicinity of the SN position using two methods. The first is 
just a simple calculation of the limiting sensitivities based 
on Poisson statistics assuming the CCD paxameters given 
above. For a 2 pixel aperture, these are V555 (3cr) = 25.3 and 
V655(5(j) — 24.8. Correcting for the CTE and applying an 
aperture correction would reduce these to V555(3a) = 25.0 
and V555(5a) = 24.5. To validate this result, we performed 
simulations of a synthetic star close to the supernovae posi- 
tion, (at a point with very similar sky values to those sur- 
rounding the possible star B). Stellar photometry using PSF 
fitting procedures within the IRAF package DAOPHOT was 
carried out, and all the stellar sources within the region of 
the supernova were subtracted from the image. The Tiny 



Tim package (Krist & Hook 1999) was used to construct a 
suitable PSF for this fitting procedure. Model PSFs of vary- 
ing magnitudes between 24.4—25.6 were made, again using 
Tiny Tim, and these were added to the "flat" image. A syn- 
thetic star at the magnitude of star B (Vsss ~25.0) is not re- 
covered with convincing significance. Placing am aperture at 
the stellar position gives a measurement of V555 ~ 25.7±0.9. 
At the theoretical 5a limit of V555 ~24.5, the object can be 
visually identified and a magnitude of 24.8±0.5 is deter- 
mined. 

We recognise that there is some statistical validity in 
the 3a detection of star B by Van Dyk et al. (2003a) and 
ourselves. However we cannot confirm, with any certainty, 
an unambiguous detection of a point source. It looks more 
likely to us that this object is a fluctuation of a fairly com- 
plex non-stellar background, which could loosely be inter- 
preted as a point source, and that the progenitor is below 
(and possibly just below) the sensitivity limit. In the rest of 
this paper we assume that we do not have a secure identifi- 
cation of the progenitor, and base our results on the robust 
5a upper limit to its V555 magnitude. We agree with Van 
Dyk et al. that there is no sign of this object in the other 
three filters, and exactly the same method for artificial star 
additions was applied to F814W and F336W, resulting in 
similar 5a detections at /8i4=24.8 and f7336=21.6. The im- 
age in the UV filter also shows nothing to a limiting magni- 
tude of mpieoBw = 19.8, and as this isn't particularly useful 
for placing constraints on the progenitor or its environment 
it is not discussed any further. 



2.2 Photometry of the stellar population around 
SN 2001du 

Photometry in the three filter images F555W, F675W and 
F814W taken after explosion on 26 November 2001 was car- 
ried out in order to help determine a reddening towards the 
surrounding stellar populations. These images were chosen 
over the prediscovery frames as they are deeper and have 
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Figure 2. A close up of the original F555W and F814W prc-explosion images before the transformation described in Section 2.1. The 
position of the SN is marked at position (0,0) with the orthogonal lines. The width of these lines represent the combined error in its 
position and the transformation RMS fit errors (0.14"). We keep the same labeling as Van Dyk et al. (2003a) for stars A-B-C, and extend 
it to the other objects detected in the field. The possible object (star B) detected at 3o" above the background noise is labeled, and its 
position is marked with a solid dot. The size of this dot represents the probable centroid error on the position of the source, and within 
the errors it is not coincident with the supernova. Star B is difficult to identify visually, and is more likely a variation in the complex 
background rather than a true point source. It is not detected in the F814W filter. The size of the image panels are 5", corresponding 
to a distance of of 460 pc at the distance of NGC1365. The directions of north and east are marked. 



better spatial resolution (as SN 2001du is on the PCI chip). 
Photometry using PSF fitting techniques was carried out 
on the stars on the PCI chip only, as the reddening to- 
wards other parts of the galaxy contained in the rest of the 
WFPC2 chips is not applicable to our present study. The 
on-the-fly recalibrated images were taken from the archive 
and the IRAF version of DAOPHOT was used to do the pho- 
tometry. Aperture photometry was first carried out with a 
4.4 pixel radial aperture. Model PSFs were calculated us- 
ing the Tiny Tim package (Krist & Hook 1999) in each fil- 
ter, and these were used in the DAOPHOT task ALLSTAR. 
The zero-points from the WFPC2 headers were applied, and 
several corrections where included. Firstly the charge trans- 
fer efficiency (CTE) correction was determined (Whitmore 
et al. 1999), then the magnitudes were corrected for the 
4.4 pixel aperture using the encircled energy tables of Holtz- 
mann et al. (1995a). Finally a colour transformation was cal- 
culated to determine magnitudes of the stars in the Johnson 
VRI bands from the WFPC2 flight system (Holtzman et al. 
1995b). Results from the photometry are discussed below in 
Section 4.1.3. 

2.3 Spectroscopy of SN 2001du 

Optical spectra of SN 2001du were taken on two epochs at 
the South African Astronomical Observatory (SAAO), and 
at the European Southern Observatory. The first spectrum 
was taken on 2001 September 2.05 at the SAAO 1.9m using 
the Cassegrain Grating Spectrograph, with the No. 6 grat- 



ing (600 lmm~^) and SITe detector providing a coverage of 
3500 — 5400A. The spectral dispersion was l.lApix"^, and 
the resolution measured from the width of the arclines was 
approximately 3A. Our second spectrum was taken at ESO- 
La Silla on the Danish 1.54m telescope on 2001 September 
9.35. The DFOSC spectrometer was employed with grism 
No. 4 and the 4kx2k EEV CCD which provided a dispersion 
of 2.8 A pix , and useful coverage between 3300 - 9200A. 
The spectra were reduced using standard techniques within 
IRAF. The bias level was subtracted from the CCD im- 
ages, before flat-fielding, extraction and wavelength calibra- 
tion from arc line maps. A photometric standard star was 
observed on each occasion for relative flux calibration. Both 
of the reduced spectra are presented in Fig. 3. 



Smartt et al. (2001c) and Wang et al. (2001) classified 

this as a typo II-P supernova from the appearance of the P- 
cygni like H i Balmer lines, particularly H/3. The unpublished 
lightcurvo of Suntzeff & Krisciunas as shown in Van Dyk et 
al. (2003b) shows it to have had a plateau phase of approxi- 
mately 100 days. The spectrum and lightcurve is very similar 
to the well studied SN 1999em which was a classic II-P event 
with a 95 day plateau (Hamuy et al. 2001; Leonard et al. 
2002b; Elmhamdi et al. 2003), and further discussion of the 
similarities are discussed below in Section 4.1.1. 
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2000 3000 4000 5000 6000 3.5 4.5 5.5 6.5 7.5 8.5 9.5 

Wavelength Wavelength / 10^ 

Figure 3. (a): The dereddened early time spectrum of SN 2001du taken on Sept. 2.05 is compared with black-body continuum fluxes 
of various temperatures (dotted lines). Comparisons with an unphysically large temperature of Tbb = lO^K effectively set a hard upper 
limit to the reddening. If the spectrum is dereddened by more than E{B — V) >0.6, then the continuum cannot be matched by any 
temperature. For comparison the more likely temperatures of the photosphere of Tbb = 12000K and 20000K (see Baron ct al. 2000), are 
shown with the reddening required to match, (b): Comparisons of the optical spectra of SN 2001du (heavy line) with SN 1999cm (lighter 
line — from Hamuy ct al. 2001). The upper spectra are SN 2001du taken on Sept. 2.05 and SN 1999em taken +7 days after explosion. 
The continuum of SN 2001du is very blue similar to SN 1999em, and indicates the reddening of the two cannot be very different. The 
difference in the slopes is likely due to the difference of around 4 days in the epochs. The lower spectra are SN 1999em taken on day +18 
after explosion, and SN 2001du taken on 2001 September 9.35 UT. The spectra are almost identical (after a linear scaling), indicating 
that SN 2001du is most probably also around +18 days, and the reddening must be similar to that for SN 1999em, for which a value of 
E{B - y) = 0.1 ± 0.05 was measured. 



3 THEORETICAL STELLAR EVOLUTIONARY 
MODELS 

We have calculated a number of stellar evolutionary models 
to compare with the luminosity constraints from the pro- 
discovery images. Models were constructed at every integer 
mass between 9— 2OM0 which is the most interesting re- 
gion probed by the observational data, plus several higher 
masses at 25, 40 and 60Mq. As described in Smartt et al. 
(2002a), the models come from the most recent version of the 
Eggleton evolution program (Eggleton 1971; Egglcton 1972; 
Eggleton 1973). The equation of state, which includes molec- 
ular hydrogen, pressure ionization and coulomb interactions, 
is discussed by Pols et al. (1995). The initial composition is 
taken to be uniform with a hydrogen abundance X = 0.7, 
helium Y = 0.28 and metals Z — 0.02 with the motooritic 
mixture determined by Anders & Grevesse (1989). Hydro- 
gen burning is allowed by the pp chain and the CNO cycles. 
Helium burning is explicitly included in the triple a reac- 
tions and reactions with ^^C, ^''N and ^^O along with carbon 
burning via ^^C +^^ C only and the disintegration of ^"Ne. 
Other isotopes and reactions are not explicitly followed. Re- 
action rates are taken from Caughlan and Fowler (1988). 



Opacity tables are those calculated by Iglesias, Rogers and 
Wilson (1992) and Alexander and Ferguson (1994). An Ed- 
dington approximation (WooUey and Stibbs 1953) is used 
for the surface boundary conditions at an optical depth of 
T = 2/3. This means that low-temperature atmospheres, in 
which convection extends out as far as r w 0.01 (Baxaffe et 
al. 1995), axe not modelled perfectly. However the effect of 
this approximation on observable quantities is not signifi- 
cant in this work (sec for example Kroupa and Tout 1997). 
We included convective overshooting in this run of mod- 
els as described in Pols et al. (1997) in order to compare 
consistently with the models of the Geneva group which in- 
corporate this mixing process. There is no mass-loss added 
to the 9 — 2OM0 tracks as this plays only a very minor role 
in the evolution of these stars. Mass loss for the higher mass 
stars is included as described in Dray et al. (2003) according 
to their NL prescription. 

The most interesting luminosity range which is of rel- 
evance to this paper, and the studies of the progenitors of 
SNe II-P, corresponds to the 9 — 2OM0 mass range. For these 
masses the end points of the evolutionary tracks at the end 
of core-carbon burning in our tracks are very close, but not 
identical, to those of the Geneva group (Schaller et al. 1992; 
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Meynet et al. 1994). The Geneva tracks are calculated at 9, 
12, 15 and 20Mo , and we have reproduced these and all in- 
teger values in between. Our tracks between 9-15Mo end at 
a luminosity 0.1— 0.15 dex higher than those of the Geneva 
group, and at 2OM0 they arc virtually identical. One should 
note that when convective overshooting is included, then 2nd 
dredge-up does not occur. Smartt et al. (2002a) discussed 
the effects of 2nd-dredge in T-llM© models without con- 
vective overshooting, which predicts that stars of this mass 
can achieve significantly higher theoretical luminosities. At 
this point however this serves as an illustrative example of 
the uncertainties in our understanding of the convective and 
mixing processes in massive stars, and the final luminosities 
of the 7-IIM0 stars cannot yet obscrvatioually constrain 
the models. This area requires more careful scrutiny from 
both the observational and theoretical side, and we intend 
to study the latter in future publications using our stellar 
evolutionary code. Given the fact that the pro-supernova lu- 
minosities in the 9-15M0 range arc approximately 0.15 dex 
higher than the Geneva tracks, in this paper we adopt the 
conservative policy of incorporating the errors on the lumi- 
nosity limits in the estimations of the upper mass limits as 
described further below. 



4 DISCUSSION OF RESULTS FOR SN 2001DU 

4.1 Reddening and metallicity of the SN 2001du 

region 

Several methods can be used to estimate the reddening 
towards the supernova, and hence towards the progenitor 
object. Those have been discussed previously in Smartt et 
al. (2001a; 2002a) and in Leonard ct al. (2002a), and here 
we try three different methods to determine the extinction. 
The method of using the Na iD interstellar absorption lines 
as a tracer of reddening along the line of sight is also of- 
ten employed for nearby SNe. However these lines are not 
visible in the Sept. 9.35 spectrum. Wo have measured the 
upper limit to their strength of EW < 0.4A which is al- 
ready suggestive of low interstellar absorption in the range 
0.06 < E{B — V)< 0.20 as deduced from the empirical rela- 
tions between the EW of the Na iD features and line of sight 
reddening given in Turatto et al. (2003). 



4.. 1.1 Comparison of the continuum slope with 
SN 1999em and a black-body function 

As discussed in Leonard et al. (2002a) and Eastman et 
al. (1996), a hard upper limit to the extinction toward a 
SN II-P in the photospheric stage can be determined. An 
arbitrarily, and unphysically, hot black-body temperature 
is compared with the dereddened continuum of the super- 
nova. If SN 2001du is doroddoned by E{B - V) >0.6 then 
even a function with Tbb = lO'^ K is unable to reproduce 
the slope of the continuum (see Fig. 3). As pointed out by 
Leonard et al. (2002a) it is desirable to have a spectrum 
taken as early as possible after explosion, because as the 
photosphere cools the upper limit that can bo derived will 
become less restrictive. Using this method they derived a 
limit of E{B - V) < 0.45 for SN 1999gi with a spectrum 
taken 0.62 days after discovery. As we discuss below, our 
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Figure 4. Two-colour diagram showing of the supergiants of 
NGC 1365, in the vicinity of SN 2001du. Also shown is the se- 
quence of theoretical supcrgiant colours (BcsscU 1990) reddened 
by E(V — I) = 0.155 ± 0.028. The reddening vector is indicated 
by the arrow. 



earliest spectrum of SN 2001du was taken approximately 8 
days after discovery. This establishes a very secure upper 
limit to the extinction, however in all likelihood the actual 
reddening will be lower, and a more applicable Tbb will be 
of the order 10000-20000 K. In Fig. 3 we plot the black- 
body fluxes along with suitably dereddened observed con- 
tinua, and conclude that the reddening towards SN 2001du 
is likely to be of the order 0.0—0.22, with a hard upper limit 
of E{B - V)=Q.6. 

The supernova SN 1999em is a very well sampled nearby 
type II-P event and has been studied in detail by Baron et 
al. (2000), Hamuy et al. (2001), Leonard et al. (2002b) and 
Elmhamdi et al. (2003). The latter three studios have exten- 
sive optical spectroscopy during the plateau phase. Hamuy 
et al. have kindly made their data available in electronic 
format on the suspect^ archive website, and wo have ac- 
cess to the Elmhamdi ct al. electronic data. The slopes of 
the contirma in the Elmhamdi et al. and the Hamuy et al. 
spectra around day -|-8 after explosion are very similar, and 
in Fig. 3 the similarity with SN 2001du is illustrated. The 
time of explosion of SN 1999om is taken as to = 2451478.8 
JD from the EPM analysis (Hamuy et al. 2001), hence the 
1999 November 14 spectrum corresponds to day -|-18. Fig. 3 
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Table 2. Limits on the bolometric magnitude and luminosity of the progenitor star of SN 2001du 
using the 5cr limits from the Vb66 and /814 images. 
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—9.08 
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— 1.78 
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—9.20 
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B8 
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AO 
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5.30 
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0.12 
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0.69 
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0.78 
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0.73 


-7.47 


1.88 


KO 
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-0.50 


-8.00 


5.10 


0.74 


-7.48 


4.89 


K2 


4310 


-0.61 


-8.11 


5.14 


0.79 


-7.41 


4.86 


K5 


3990 


-1.01 


-8.51 


5.30 
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-7.11 
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MO 
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4.85 
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-3.47 


-10.97 


6.28 


0.67 


-7.53 


4.91 



shows the excellent match between the +18 day spectrum of 
SN 1999em and our SN 2001du ESO spectrum taken on 2001 
September 9.35. We attempted to match the SN 2001du 
spectra with the Hamuy ot al. data for +13 (November 9) 
and +23 (November 19) days and applied variable redden- 
ing to the spectra of SN 2001du but could not consistently 
match the continuum slope, the strength of the absorption 
lines and the Ha P-Cygni feature. In addition we varied the 
absolute linear scaling factor applied to the SN 2001du spec- 
trum (to account for distance and intrinsic magnitude dif- 
ferences) and the amount of dereddening (assuming R=3.1) 
when matching it to the +18 day SN 1999em spectrum. We 
found that any extra reddening added to SN 2001du was 
incompatible with a good match and oven applying an addi- 
tional reddening of E{B — V) = 0.1 gave poor results. Hence 
the spectrum of SN 2001du taken on 2001 September 9.35 
UT is virtually identical with that of SN 1999em at +18 days 
after explosion and implies an identical reddening. Baron et 
al. (2000) have compared the early time optical and UV 
spectra of SN 1999em with a non-LTE model atmosphere 
fit. They simultaneously determine a temperature and a low 
value of reddening of E{B — V) — 0.10 ± 0.05, and argue 
that increasing the value of E{B — V) will not lead to a con- 
sistent fit to all the features in the spectra. In summary, the 
continua and line strengths of SN 1999em and SN 2001du 
on +18 days are virtually identical, and any further artificial 
reddening of SN 2001du leads to significant inconsistencies. 
Hence our estimate of reddening towards SN 2001du using 
this method is E{B -V)=0.1± 0.05. 



4- 1.2 Estimates of reddening from nearby Hii regions 

The abundance gradient in NGC1365 has been studied by 
Roy & Walsh (1997) who observed the nebular lines of 55 H II 
regions. For each Hll region they determined an extinction 
at H/3, which is the combined value of Galactic reddening 
{Av = 0.16), and extragalactic extinction. Their nebular 
region RW21 is virtually coincident with the star formation 
region around the position of the supernova, and the cen- 
tre of this Hii region is estimated to be approximately 4" 
from SN 2001du. Hence it is likely that this nebulosity is 
directly associated with the ionizing hot stars in the vicinity 
of SN 2001du i.e. the bright clustered sources in the Usse 
image in Fig. 1(a). Roy & Walsh (1997) measured a log- 
arithmic extinction at 486lA of c(H/3) = 0.83 ± 0.06 for 
RW21 and using an extinction law of A{-i86l) / A{V) = 1.16 
(CardelU et al. 1989) this translates into Ay = 1.8 ± 0.1, 
or E{B -V) = 0.58 ± 0.03 (assuming R=3.1). There are 
two other Hii regions lying in the same spiral arm at the 
end of the central bar in NGC1365 both at distances of ap- 
proximately 15" from the supernova (RW19 and RW20). 
These have a virtually identical extinction of E{B — V) = 
0.53 ± 0.03, hence the extinction toward the ionizing regions 
near the SN position which produce the nebular lines is of or- 
der E{B-V) = 0.58 + 0.03. However as discssed by Calzetti 
et al. (1994) and Calzetti (1997) the extinction measured 
in starformiug regions from the nebular linos is generally a 
factor of two higher than that of the underlying stellar con- 
tinuum. Hence the H ii region reddening is not likely to be 
directly applicable to the SN 2001du reddening as the super- 
nova it is not coincident with the position of the UV-bright 
ionizing O-stars. We hence assume that a reddening a fac- 
tor of two lower should be applicable to the stars near the 
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RW21 region i.e. E{B -V) = 0.29 ± 0.03. While this may 
appear somewhat arbitrary, we should get a reliable mean 
when combining it with the other extinction estimates . 

Roy & Walsh determined the oxygen abundance in 
their Hll regions from an analytic calibration of the R23 
ratio ii23=([OII]A3727+[OIII]AA4959,5007)/H/J. While this 
method usually leads to a reasonably consistent estimation 
of how abundances differentially vary within a galaxy, the 
absolute abundance determined for any particular J?23 ratio 
is highly dependent on the calibration employed. Recently 
Smartt et al. (2001b) and Trundle et al. (2002) have dis- 
cussed various literature calibrations of the R23 ratio, and 
compared extragalactic nebular abundances directly with 
massive stars born within the H 11 regions. They suggest 
that the calibrations of Pilyugin (2002) and McGaugh (1991) 
(the latter parameterized by Kobulnicky et al. (1999)) give 
the best agreement for oxygen abundances in the gas and 
in massive B-type supergiants. Applying these calibrations 
to the _R23 ratio of RW21 results in an oxygen abundance 
12 + log (0/H)=8.5dex (from the Pilyugin calibration), and 
8.8 dex (from McGaugh). Hence within the uncertainty of 
the absolute calibration of this method, the H 11 region RW21 
is approximately solar in mctallicity; assuming the solar oxy- 
gen abundance is 8.83 dex (Grevesse & Sauval 1998). 



4.1.3 The reddening towards surrounding stars in 

NGC1365 

In this case the reddening can be measured from 3-colour 
photometry of the bright supergiants in the area around the 
supernova, and one then assumes that this is representa- 
tive of the reddening toward the supernova itself. Photom- 
etry in three filters (F555, F675W and F814W transformed 
to standard VRI) was measured for 190 stars on the PC 
chip in the vicinity of SN 2001du (to a radius of 23"from 
the SN position) . The locus of the observed stars, in a two- 
colour plane {V — I ,V — R) was compared with a theoretical 
supergiant colour sequence from Bessell (1990) for galactic 
supergiants. The reddening was determined by ascertaining 
the displacement between the observed stars and the the- 
oretical sequence in the two-colour plane. A weighting fac- 
tor was included in the calculation, in favour of those stars 
with more accurately determined magnitudes. The redden- 
ing was determined to be E{V - I) = 0.16 ± 0.03. Fig. 4 
shows the two-colour diagram of the stars imaged in the 
vicinity of SN 2001du; also shown is the theoretical super- 
giant sequence appropriately reddened. Utilising a redden- 
ing law of Rv = Av/E{V - /) = 2.45 (Silbermann et al. 
1999) yielded a value of Ay = 0.38. In determining a red- 
dening towards the Cepheid population of NGC1365 (in or- 
der for distance determinations), Silbermann et al. (1999) 
measured a value of Av = 0.40. Within the errors this is 
identical to our measurement in the field of SN 2001du, and 
suggests that the reddening towards this part of the galaxy 
is not significantly different to that derived by Silbermann 
et al., and is not prohibitively high. A standard reddening 
law of Rv = Av/E{B — V) = 3.1 would imply a value of 
E{B — V) = 0.12 ± 0.02 is appropriate for the Une-of-sight 
to SN 2001du from this method. 



4.2 A luminosity and mass limit for the 
progenitor of SN 2001du 

The three methods above give reassuringly consistent mea- 
surements for the reddening towards SN 2001du, and pre- 
sumably its progenitor star. Hence we adopt a mean value 
of these estimates of E{B — V)= 0. 17 ± 0.09 and (assuming 
Rv = 3.1) Ay = 0.53 ± 0.26. Smartt et al. (2001a; 2002a) 
discuss the general technique employed to determine an up- 
per luminosity limit for the progenitor for it not to have 
been detected on prc-cxplosion images. The pre-explosion 
?7336 archive image was not sufficiently deep to place any sig- 
ruficantly constraining upper luminosity limits for the pro- 
genitor star and is not considered further in the following 
discussion. The V555 and /gi4 5a detection limit magnitudes 
were transformed to V = 24.4 and / = 23.5 using the trans- 
formations of Holtzman et al. (1995b). The corresponding 
absolute magnitudes {My and Mi) were calculated by cor- 
recting for the extinction, determined earlier, and adopting 
a distance modulus to NGC 1365 of (ni-M) = 31.39 ±0.26 
(Ferrarese et al. 2000). The limiting absolute magnitudes 
were determined to be: My = —7.5 and Mi = —8.2. Bolo- 
mctric magnitudes, for each waveband, were determined by 
applying bolomctric corrections, and an appropriate colour 
correction {V — I)o for Mi, both taken from Drilling & Lan- 
dolt (2000). The luminosity limits, as a function of tem- 
perature, were calculated from the bolometric magnitudes 
assuming a solar bolometric magnitude of Afboi = +4.74, 
and are listed in Table 2. These luminosity limits were over- 
laid on an HR diagram with our evolutionary tracks (as 
described in Section 3) . An HR diagram showing the exclu- 
sion region defined by the luminosity upper limits is shown 
in Fig. 5. A closer view of the end-points of the tracks in 
the red supergiant region, and the exclusion limit defined 
by the I band sensitivity limit is further shown in Fig. 6. 
The expected magnitude errors in distance modulus (±0.26), 
reddening (±0.26), sensitivity limits (±0.2), and BC and 
intrinsic stellar colours (±0.4) are combined in quadrature 
giving an uncertainty on Mboi of ±0.58 mag. As we pointed 
out in Sections, the end-points of our tracks in the red- 
supergiant region tend to be systematically higher than the 
Geneva tracks by 0.1— 0.15 dex. As the observational uncer- 
tainties are significant we adopt the conservative approach 
of adding the total error to the luminosity upper limit to de- 
rive a robust upper mass- limit. In this case we can say that 
the upper limit to the initial mass of the progenitor star of 
SN 2001du is 15M0. 

As discussed in Sect. 2.1 we believe it unlikely that star 

B is a single point source which is precisely coincident with 
the the supernova and is an unlikely progenitor candidate. 
However if it was the progenitor it would have My = —6.8, 
and V — I < 1.4. This would constrain it to be just below 
the exclusion region in Fig. 5 and to be bluer than approx- 
imately K2 spectral type i.e. cooler than logTcff <^ 3.7, and 
with a mass in the region 15-2OM0. Hence it would not 
have been an M-type supergiant, but an unperturbed early 
K-type could provide a large enough hydrogen envelope to 
sustain a plateau pheise of ~100 days (see Table 3), and is 
thus not inconsistent with the supernova evolution. Van Dyk 
et al. (2003b) have also reanalysed the site of SN 2001du us- 
ing our images of the supernova, and come to the same con- 
clusion that Star B is not the progenitor. The upper mass 
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Figure 5. SN 2001du: Luminosity limits from pre-explosion im- 
ages in the V and / are plotted as thick solid lines. The shaded 
region is where a progenitor would have been detected in at least 
one filter. Overlaid are the evolutionary tracks as described in 
Sect. 3 for stars with main sequence masses 9 — GOMq, for so- 
lar metallicity. The locations of the progenitors of SN 1987A and 
SN 1993J are also shown. 




Figure 8. SN 1999em; updated exclusion region diagram for the 
progenitor of SN 1999em (from combined R and I limits), using 
the new distance estimate for NGC1637 of llMpc. Our evolu- 
tionary tracks as discussed in Sect. 3 are plotted along with the 
uncertainties on the luminosity limits. The new distance increases 
the upper mass limit to ISMq. 




Figure 6. SN 2001du; End stages of the evolutionary tracks for 
stars with initial masses in the range 9— 2OM0. The limiting lu- 
minosities are shown as the thick solid lines. Uncertainty limits, 
on the I-band limiting luminosity are shown as the dashed lines. 
The limiting luminosity including the expected error constrains a 
robust upper mass limit of 15Mq 

limits derived, again repeating our method, are in agree- 
ment with those derived here (13M<^4), although we have 
reduced the uncertainties using our finer grid of stellar evo- 
lutionary tracks. 



5 DISCUSSION OF RESULTS FOR TWO 
OTHER II-P SUPERNOVAE 

5.1 SN 1999em 

In Smartt et al. (2002a) we discussed the prediscovery im- 
ages of SN 1999em which wore taken from the ground with 
the Canada- France-Hawaii telescope at a resolution of 0.7". 
As with SN 2001du and SN 1999gi (see Section 5.2) there 
was no detection of a progenitor star at the position of 



SN 1999em, and the sensitivity limits of the VRI frames 
allowed an upper mass limit to be derived. Similar meth- 
ods as discussed above led us to determine an upper mass 
limit of approximately 12M0. However these conclusions are 
based on the assumption that the ground-based images do 
not have contaminating objects within the seeing PSF of 
0.7". As shown in the case of SN 1999gi, which exploded in a 
fairly large OB-association, it is quite possible that ground- 
based resolution images would not be adequate to resolve 
the individual massive stars. For that reason we have reob- 
served the environment of SN 1999em with WFPC2 dur- 
ing our Cycle 10 program GO9041. The main goal of this 
observation was to investigate the immediate environment 
surrounding the progenitor of SN 1999em, and determine if 
a contaminating star cluster would compromise the result. 
SN 1999em was observed with WFPC2 on 2002 January 15 
through the filters F555W, F675W and F814W. The F555W 
image is shown in Fig. 7, and we have identified the posi- 
tion of the supernova through matching the brightest single 
stars in the field of the WFPC2 data and the WHT im- 
age of Smartt et al. (2002a). A geometric transformation 
between the two frames results in a residual uncertainty of 
approximately 0.1" in the SN position. The WFPC2 im- 
age shows a faint object at the position of the supernova, 
which is not a clean point source. Two peaks appear, super- 
imposed on a quite faint extended region of flux. Although 
this diffuse emission is faint it is clearly visible in both the 
F555W, and the F675W but not the F814W filter. We esti- 
mate the position of the brighter peak at the centre of Fig 7a 
is within 0.08" of the supernova centroid, which is within 
the errors of our differential astrometry. The extended weak 
flux is approximately 0.5" in diameter, or 25 pc (assuming 
the distance discussed below). This is much too far to be 
due to a light echo from SN 1999cm and is likely a small 
star cluster which hosted the progenitor star. The magni- 
tude of this object as a whole (from aperture photometry 
using a radius of 0.3") is V555 = 23.7 ± 0.2, which is a com- 
bination of the residual supernova magnitude and the un- 
derlying association. This would imply that the host cluster 
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Figure 7. (a): A WFPC2 F555W image of the SN1999em taken on 2002 January 15, some +810 days after explosion. The position of 
the supernova was identified from the high-resolution WHT images presented by Smartt et al. (2002a). (b): Close up of the faint source 

at the position of SN 1999cm, which is not a simple point source (compare with the other two point sources at j/ ~ 0.6). It appears as 
two peaks superimposed on a faint diffuse background. However the region around the SN position is uncrowded and relatively clear of 
contaminating bright sources within the 0.7" seeing disk of the prc-explosion frames of Smartt et al. (2002a). This validates the use of 
these lower resolution frames in setting a luminosity and mass limit for the progenitor star. 



has Mv > —6.1, which is below the detection limit of the 
prc-cxplosion frames {Adv ~ —6.49; Smartt et al. 2002a). 
Hence the size and luminosity of the cluster is fairly typical 
for the smallest young associations found in spirals and in 
the Local Group e.g. BresoUn et al. (1996; 1998). For this 
paper, the most important conclusion is that the images con- 
firm that within the seeing disk of our prediseovery CFHT 
images there are no other bright contaminating objects. 

Leonard et al. (2003, in preparation) have observed the 
galaxy NGC1637 on several epochs with HST/WFPC2 dur- 
ing Cycle 10 to measure a Cepheid distance to the galaxy. 
This is the first attempt to compare a Cepheid distance 
to the expanding photosphere method as applied to the 
type II supernovae, and the Cepheid distance measurement 
is D ~llMpc (D. Leonard and N. Tanvir, private com- 
munication). Hence we have revised the exclusion diagram 
presented in Smartt et al. (2002a) with this new distance. 
As the HST/WFPC2 images confirm that no contaminat- 
ing stellar sources exist within the ground-based seeing disk 
that would compromise the point source sensitivity limits, 
the diagram shown in Fig. 8 should now be a robust result. 
Applying the error estimates and (as in Smartt et al. 2002a) 
the same conservative criteria as applied to SN 2001du of 
using the highest mass that falls within the upper bounds 
of the errors, we now derive an upper mass limit of ISM©. 
The fact that the underlying association or cluster, in which 
the progenitor was born, was small adds further weight to 
our argument that the progenitor could not have been very 
massive. Assuming a typical initial mass function, it is sta- 
tistically unlikely that a cluster of this size would have given 
birth to very massive stars of greater than 2OM0 . The new 
limit derived is tabulated in Table 3 for comparison with 
other known progenitors. 




Figure 9. SN 1999gi: End points of stars with initial masses in 
the range 8— 2OM0 using our evolutionary tracks as described 
in Sect. 3. The limiting luminosities (from the F606W filter) are 
shown as the thick solid lines and the uncertainties on these are 
the dashed lines. This gives a robust upper limit of 12M0. 



The unexpected faint magnitude of the supernova has 
implications for the Elmhamdi et al. (2003) conclusions. 
Even assuming the total magnitude measured (i.e. super- 
novae plus underlying cluster) to be an upper limit to 
the SN 1999em flux itself, it is significantly fainter than 
the last reported detection by Elmhamdi et al. (2003) of 
V = 20.05 ± 0.18 on 2001 March 16. Elmhamdi et al. 
have suggested that their last photometric point is fainter 
than the expected linear decay (extrapolated from previous 
points) of 7„ ~ 0.97 mag(100d)~^ by approximately 30%, 
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Figure 10. Estimate of the lower mass limit for the progenitor of 
SN 1997bs. The object was detected at My =i -8.1 (in F606W; 
Van Dyk et al. 2000), and the bolometric luminosity locus as a 
function of Tefj (from the supergiant table in Drilling and Landolt 
2000) is plotted. Overlaid are evolutionary tracks for stars with 
main sequence masses 12 — 40Mq, for solar mctallicity. The pro- 
genitor must have been somewhere on this line prior to explosion, 
hence an approximate lower limit to the mass is M> 20Mq. 



and have interpreted this as a sign of dust formation. If 
we extrapolate the linear decay of Elmhamdi et al. (from 
the last point on 2001 February 01), then the SN 1999em 
should be V~23.0 (at day +810) in the HST WFPC2 im- 
ages. As it is O.?"* fainter (assuming the F555W bandpass 
is close to V) this is supporting evidence that the devia- 
tion they reported is real, and continues even more severely 
at late times. Elmhamdi et al. also suggested that the [O i] 
line profile variations after day ~465, were indicative of dust 
formation in the cooling supernova. This faint WFPC2 mag- 
nitude of the SN 1999em would seem to add further support 
to that. The multi-epoch images of Leonard et al. (GO 9155) 
to monitor the Cepheid variable population also contain the 
SN 1999em on the field of view, hence this will allow them 
to probe the late time lightcurve in F555W and F814W, and 
dust formation scenario, in much more detail. 



5.2 SN 1999gi 

In Smartt et al. (2001a) we presented the HST archive pre- 
explosion images for SN 1999gi along with images of the 
supernova taken 14 months later in order to astrometrically 
determine the supernova position precisely. There was no de- 
tection of a progenitor star in either the F300W or F606W 
filter, and using the methods described above in Section 4.2 
we derived an upper mass limit for the progenitor star of 
9^2 M©, which assumed it weis a red supergiant in the final 
stages of its evolution. However in this paper we assumed a 
distance of 7.9 ± 2 Mpc, and recently Leonard et al. (2002a) 
have determined an improved distance to SN 1999gi using 
the expanding photosphere method. They estimate a dis- 
tance of 11.1 ± 2 Mpc, and have also carried out a careful 
comparison of reddening towards the supernova from four 
different methods (which we followed in Sect. 4.1). Leonard 
et al. estimate a slightly higher value of E{B — V) = 0.21 
compared to the E{B — V) = 0.15 used in Smartt et al. 



(2001a), and derive an upper mass limit for the star of 
I5I3MQ on this basis. 

We have recalculated the luminosity and mass limits 
using the new distance of Leonard et al., and the higher red- 
dening, but using the better sampled tracks at our disposal 
we can reduce the uncertainty in the mass range of Leonard 
et al.'s result. Leonard et al. have estimated the mass limit 
by adding the upper error uncertainty and simply taking 
the nearest modeled progenitor mass to that, with an upper 
limit of the next highest mass model in the Geneva tracks. 
The exclusion region of the HRD is plotted in Fig. 9, and 
applying the same criteria as for SNe 2001du and 1999em 
we get an upper mass limit of I2M0, which is below our 
best estimate and the uncertainty. The difference between 
this limit and the ISIjMq estimate by Leonard et al. is 
that our tracks end at slightly higher luminosities than the 
Geneva models, and also we can rule out the rather large 
range by employing the integer sampled mass tracks. It is 
understandable that Leonard et al. employed a conservative 
approach in estimating the upper limit and errors. However 
we can certainly rule out the progenitor being a 2OM0 red 
supergiant, as this lies significantly above the detection lim- 
its it would have a luminosity of log L/L© ~ 5.3, and 
hence a magnitude of Vsoe — 23.9, which would have been a 
clear (lOcr) detection on the pre-explosion images. In addi- 
tion there may be a slight discrepancy between our temper- 
atures of pre-supernova M-supergiants and those in Fig. 10 
of Leonard et al., where the effective temperatures of the 
evolutionary track endpoints appear too cool by ~0.05 dex. 
They use the Geneva tracks (as employed in Smartt et al. 
2001a), but the plot does not exactly match the tabulated 
values i.e. their end-points would be close to logTo£F=3.5 in 
our Fig 9. If the tracks are shifted to the tabulated temper- 
ature then their luminosity limit of the F606W filter would 
no longer skim the 2OM0 track, and the 15Mq track would 
be inside the exclusion region. 

We accept the distance and reddening determinations of 
Leonard et al. as more appropriate to those used in Smartt 
et al. (2001a), and that Leonard et al. are being understand- 
ably cautious in using the estimate plus uncertainty to de- 
rive a very hard upper limit. But the simplest interpreta- 
tion of Fig. 9, where we use better sampled end points in 
the tracks, indicates an upper mass-limit of 12M0. Clearly 
these Eissumptions rely on how accurate the theoretical stel- 
lar evolutionary models are in predicting the luminosity and 
temperatures of M-type supergiants. 

5.3 Progenitors of core-collapse supernovae 

In this section we review all of the direct observational evi- 
dence that exists for core-collapse precursors. The events for 
which we have some restrictive information on the progeni- 
tors are listed in Table 3. Only those supernovae which are 
fairly certain to be core-collapse explosions, and have a well 
defined SNe type attributed are included. 

The event SN 1978K is not included although it had a 
progenitor object visible. The supernova itself was serendip- 
itously discovered in January 1990 as a peculiar object on an 
archive photographic plate (Ryder et al. 1993). Hence there 
is only a crude light-curve from archive plates taken cover- 
ing that region of sky. Ryder et al. (1993) conclude that it 
was more likely to be a type II event rather than an LBV 
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Table 3. A list of known supernova progenitors and constraints on their masses and spectral types for supernovas which are certain 
to be real core-collapse events and have been reliably typed. The estimated initial mass of the progenitor star is listed as Mj and is 
from direct measurements of the progenitor star or limits thereon. Compilation values for mass and metallicity are taJcen from Smartt 
et al. (2002a), and updated with results from this paper and Smartt et al. (2002b; for SN 2002ap). We estimated the metallicity of the 

progenitor of 1997bs simply from its distance from the centre of NGC3627, and applied a typical abundance gradient for an Sb galaxy 
from Vila-Costas & Edmunds (1992). The mass of the ejecta in the supernova is Mej and the radii of the progenitors are estimated from 
analysis of the supernova spectra and lightcurve, as described in Section 5.3. The values for SN 1987A and SN 1993J, and SN 1980K are 
taken from Arnett (1996) and Woosley et al. (1994). 
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Figure 11. The theoretical SNe populations from the models of Heger et al. (2003), but ignoring their extremely metal poor and 
metal-free objects. All objects below the dotted line should be weak SN due to black-hole formation by fall-back of most of the material 

surrounding the iron core. All objects to the right of the solid black line should have lost their hydrogen envelope at the point of 
core-collapse, and hence should be hydrogen free SNe. The blank white regions are where no supernova is seen; below 9M0 due to the 
formation of white dwarfs in low mass stars; above 4OM0 due to direct black hole formation. The observed positions of the supernovae 
progenitors in Table 3 are approximately estimated and placed, with arrows indicating those with limits. To place the stars in the 
metallicity dimension we have assumed a zero-point on the y-axis of log Z/Z0=-1. 
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"super-outburst", however this is not beyond doubt. If it 
was a supernova we do not know its type, peak magnitude 
or very much about the hghtcurvo. Given the limited reso- 
lution of the photographic surveys (particularly close to the 
plate limit) wo do not know if the object at the position of 
SN 1978K is really a point source, a nebula or cluster (Ry- 
der et al. quote a distance of 4.5Mpc). These uncertainties 
lead us to omit SN 1978K from Table 3 for now, at least 
until some of the questions can be definitively answered. 
The event SN 1961V is also not included for similar rea- 
sons. It is very unclear if it was a genuine supernova or an 
S-Doradus outburst from a very massive star like 77 Carinae 
(Humphreys et al. 1999). The progenitor object was possibly 
identified as a star with an absolute magnitude Mboi — ^11, 
which is much too bright for a stable massive star or 77 Car 
like object in a quiescent state. Hence it was likely to be in 
outburst during its period as an Mboi — —11 object, and 
given the sparse data available we do not have a definitive 
model for the physics behind the eruption. Again because of 
these uncertainties, we do not discuss this further. 

A detection of a progenitor star associated with 
SN 1997bs in M66 has been claimed by Van Dyk et al. 
(2000), who originally suggested that this event was also 
not a genuine core-collapse supernovae but another super- 
outburst similar to SN 1961V and 77 Car, and that the pro- 
genitor actually survived the explosion in HST images. How- 
ever Li et al. (2002) have observed SN 1997bs at a later 
epoch with HST and observed it to fade considerably, and 
get bluer, compared to the Van Dyk et al. magnitudes. With 
this new information, the conclusion is that SN 1997bs was 
indeed likely to have been a subluminous Iln supernova, and 
the progenitor has not survived. Van Dyk et al. (2000) have 
estimated the absolute magnitude of the progenitor star to 
be Mv ~ -8.1 from the single exposure F606W WFPC2 
image. As no colour information is available we do not know 
the temperature or appropriate bolomotric correction to al- 
low the progenitor to be approximately placed on an HR 
diagram. However given this detection in a single band, a 
lower limit can be placed on the mass of the progenitor. In 
Fig 10 we plot the locus of a My = —8.1 star as a function of 
effective temperature and hence bolometric correction. The 
progenitor must have existed at some point on this locus, 
hence allowing a lower limit to its mass of approximately 
2OM0. 

The brightest SNe Ic to occur in recent times was 
SN 2002ap in M74 (=NGC628), and its early lightcurve and 
spectral evolution was recently studied by Mazzali et al. 
(2002). By comparing explosion models of CO cores com- 
bined with spectral synthesis to model the optical obser- 
vational data, they estimated an ejecta mass of between 
2.5— 5M0. With an assumed (but uncertain within a fac- 
tor ~2) remnant mass of roughly 2.5Mo, this would imply 
the mass of the CO precursor star should be in the con- 
servative range 2.5— lOM©. The final masses of Wolf-Rayet 
stars of the WC type in the LMC have been estimated by 
Crowther et al. (2002) to be in the range 11 — I9M0, and 
those in the galaxy in the range 7— 14M0 . Given the uncer- 
tainties, there is reasonable agreement between these two 
estimates suggesting a WC type progenitor (of initial mass 
^SOMq) is consistent with the characteristics of SN 2002ap. 
Deep ground-based images of the pre-explosiou site of the 
supernova were presented by Smartt et al. (2002b), and no 



detection of a massive star was found. This paper ruled out 
most of the upper HR-diagram as possible progenitor sites, 
although this is to be expected as SNe Ic show no evidence of 
the jjrogcnitors having H-rich atmospheres. It is likely that 
the precursor was a WR star in the WC phase, and had 
an initial mass of less than ^40Mq . However initially lower 
mass stars which have had their H-envelopes stripped dur- 
ing mass-transfer in a binaxy system, rather than through 
a radiatively driven stellar wind, still do remain viable pre- 
cursors. 

The three II-P supernovae discussed above (SNe 
2001du, 1999em, 1999gi) are also hsted in Table 3, along 
with a similarly determined upper mass limit for SN 1980K. 
Pre-explosion images for this supernova were presented by 
Thompson (1982), and we have recalculated the limit using 
the same evolutionary tracks as employed above. The more 
definite values for SN 1987A and 1993J are included. The 
analytical model for the plateau stage of type II supernovae 
of Popov (1993) produces equations relating the observed 
quantities of tp (the plateau duration), My (absolute visual 
magnitude at a representative point during the plateau), and 
Mphot (the expansion velocity of the photosphere) to three in- 
trinsic properties of the supernova. These are energy of the 
explosion and the mass and radius of the initial envelope 
of the progenitor star. More detailed hydrodynamic mod- 
els of SNe II-P by Litvinova & Nadyozhin (1985) have pro- 
duced similar relations, but with different cocfBcients. The 
mass and initial radius of the ejected stellar envelope for 
1999em and 1999gi has been estimated by Elmhamdi et al 
(2003), Hamuy (2003) and Nadyozhin (2003) using both sets 
of equations but differing values for tp, My and Uphot- The 
mass and radii derived in each paper differ quite substan- 
tially, depending on what values are assumed. The main dif- 
ferences are the distance assumption, which critically affects 
My, and the epoch at which Mphot is measured. We have re- 
calculated the range in feasible values using the most up to 
date values for distance and reddening (as discussed above), 
and get values of My = -16.7 (1999em) My = -15.9 
(1999gi) at approximately 50 days after discovery (i.e. in the 
middle of the plateau). We assume that tp ~ 100 days and 
vary Uphot between 3000— 4000km s~^ . There appears no 
discrepancy in the actual measurements of Mphot in 1999em 
between Elmhamdi et al. (2003), Hamuy et al. (2001) and 
Leonard et al. (2002b), but there is some difference as to 
which epoch is chosen to measure Mphot to use in the cal- 
culations. The ranges of ejected mass (Mej ) and progenitor 
radius (R) are tabulated in Table 3. If we allow '^2Mo to be 
lost by mass-loss and neutron star formation, then the mass 
ranges for the envelope are in reasonable agreement with 
our direct estimates for the initial mass of the progenitor 
stars. The large range in radius of the envelope is consis- 
tent with spectral types between GO— MO. We do not have 
enough photometric information available for SN 2001du to 
estimate those values for comparison. However our direct 
mass limit for SN 2001du is similar to those for SNe 1999em 
and 1999gi. Given the very similar spectra of SN 2001du 
with both of the latter (see Fig. 3), this is not altogether 
surprising. This supports our previous suggestion (Smartt et 
al. 2002a) that normal SNe II-P result from fairly low mass 
stars (<12M0) with large radii in the M-supergiant phase 
that have undergone only moderate mass-loss. The fact that 
1999em and 1999gi were faint in the x-ray and radio suggests 
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that their progenitor stars had undergone fairly moderate or 
low mass-loss. The fluxes in these wavebands have allowed 
Schlegel (2001) and Pooley et al. (2002) to estimate mass- 
loss rates for the progenitors of these two supernova which 
arc the order of M --^ 1O~^M0 yr^"*^, with a conservative un- 
certainty of a factor of two. This supports our result that 
the progenitors were not very massive stars, as one would 
expect such stars to have gone through more drastic mass- 
loss events and hence to have produced higher x-ray and 
radio fluxes. Hence wo conclude that moderate mass stars 
(of less than approximately 15Mq) which end their lives in 
the red supergiant phase tend to give rise to a fairly homo- 
geneous class of SNo II-P , and initially higher mass stars 
likely give rise to more heterogeneous supernovae, although 
clearly more results are needed to increase the statistics. 

Recent theoretical work by Heger et al. (2003a; 2003b) 
on the the pre-supernova evolution of massive stars as a 
function of mass and mctallicity, has produced supernova 
population diagrams. These relate the initial mass of a star 
to the type of supernova that is produced, and how that 
roughly depends on mctallicity of the original star. The 
metallicity plays a key role in determining the mass-loss rate, 
which in turn is a critical factor in the evolution. This has 
initially been done for single stars, with no rotation consid- 
ered. We have reproduced the figure of Heger et al. (2003b), 
and ignored the very metal poor and metal-free regions as 
they are not applicable to any of the local Universe super- 
novae discussed above. The regions between the H-P, H-L 
and Ib/c regions are delineated by the mass of the hydro- 
gen envelope which is left just before core-collapse, which in 
turn is dependent on mass-loss history and hence mctallic- 
ity. In Fig. 11 we have put the positions of the supernovae 
listed in Table 3 to compare with the model predictions. We 
emphasise that at present this is a rather qualitative com- 
parison given the approximate metallicity scale proposed on 
the ordinate, and the limits wo have for the various nearby 
supernovae. Nevertheless it is certainly valid and essential 
to compare the state-of-the-art observational and theoret- 
ical results in this data starved area. The positions of the 
three SNe H-P are in good agreement with the models, as 
is the position of SN 1987A. However SN 1980K was a type 
II-L but appears in the theoretical H-P region. Even if the 
metallicity is uncertain by a factor of two, it would be dif- 
ficult to envisage pushing this into the II-L region as it 
stands. This suggests that the progenitor lost more mass 
than one would expect for a 2OM0 star of metallicity be- 
tween 0.5— IZq. The extra mass-loss could be accounted for 
by interaction within a binary system, with the progenitor 
losing its mass through mass-transfer rather than radiatively 
driven winds. SN 1997bs is also likely to have been too mas- 
sive for the II-P region, although given that we have only a 
iower limit, the actual mass could sit comfortably in the II-L 
area. One serious discrepancy is obvious in that SN 2002ap 
was a bright Ic supernova, however the models suggest that 
these only come from either metal-rich stars, or very massive 
stars with Mi > 6OM0 . As discussed in Smartt et al. (2002b) 
the initial mass is unlikely to have boon as high as ~6OM0, 
hence there is a discrepancy with the position of the weak 
SN region. Smartt et al. (2002b) estimated the mctallicity to 
be roughly O.5Z0 , and again even if this was solar the posi- 
tion of SN 2002ap would still be discrepant with the normal 
SN Ib/c population region. Prom the Heger et al. figure one 



would expect that bright Type Ib/c only come from massive 
single stars if the metallicity is quite high, and hence they 
should be relatively rare. However invoking mass-loss again 
through close binary evolution at low metallicities could ac- 
count for objects such as SN 2002ap. The future statistics 
of masses and mass limits for the progenitors of SNe Ib/c 
will be of crucial importance in distinguishing between the 
two scenarios. In summary then it would appear that the 
theoretical II-P region stretches to masses which are a little 
too liigli. and that the weak SN region should be confined to 
lower metallicities or much higher masses. This is of course 
a very preliminary comparison, with the caveat that only 
single stellar evolution is included in the supernova popu- 
lation diagram. The prospects of setting more constraining 
limits on SNe progenitors or detecting the stars themselves 
will allow better observational constraints to be placed on 
these models in the future (see Smartt et al. 2002a for a 
discussion of future prospects) . 

6 CONCLUSIONS 

We have investigated HST images of the site of the neaxby 
supernova SN 2001du taken 6.6 years before explosion, and 

compared our results with similar pro-explosion investiga- 
tions of other core-collapse supernovae. Our findings are 
summarised below 

(i) Although there is a 3a detection of a source very close 
to the supernova position in the WFPC2 V^— band image, 
it is not precisely coincident with the supernova centroid. 
It is slightly beyond the error radius in our analysis and 
is not a clear detection of an unambiguous point source. 
We conclude that the progenitor star was likely below the 
detection limit of the pre-discovery images and is not visible 
on any of the frames. Using the 5a sensitivity limits of the 
images we estimate that the progenitor likely had an initial 
mass of <15M0. 

(ii) SN 2001du was a type H-P, and probably very similar 
to two other nearby recent SNe II-P 1999gi and 1999em. We 
have previously set similar luminosity and mass limits on 
the progenitors of those events (Smartt et al. 2001a; Smartt 
et al. 2002a) and in this paper we have revisited those us- 
ing new data. A new distance to NGC3184 determined by 
Leonard et al. (2002a) slightly changes the mass limit for 
SN 1999gi to approximately 12M0. New HST images of the 
site of SN 1999em confirm the validity of our previous re- 
sults, and using another new distance increases this upper 
mass to approximately 15M0. 

(iii) Hence all three SNo had very similar mass progenitor 
stars, which is not a surprising result given the photometric 
and spectroscopic similarity of the supernovae themselves. 
The results are consistent (but not uniquely definitive) with 
the idea of SNe II-P arising in red supergiant progenitors of 
moderate mass. 

(iv) We have made first attempts at comparing all the 
known, direct, information available on supernova progen- 
itors with theoretical models of stellar evolution up to the 
core-collapse stage. These models have produced supernova 
population diagrams which relate the initial mass and mctal- 
licity of massive stars to the types of supernovae that they 
produce. We find encouraging reasonable agreement between 
the observational and theoretical results for the masses of 
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II-P supernovae progenitors. However we find discrepancies 
for a II-L and briglit Ic event. This is preliminary work that 
should increase in statistical significance over the coming 
years 
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